Dopant Spatial Distributions: Sample Independent Response Function And Maximum 

Entropy Reconstruction 
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We demonstrate the use of maximum entropy based deconvolution to reconstruct boron spatial 
distribution from the secondary ion mass spectrometry (SIMS) depth profiles on a system of vari- 
ously spaced boron 5-layers grown in silicon. Sample independent response functions are obtained 
using a new method which reduces the danger of incorporating real sample behaviour in the re- 
sponse. Although the original profiles of different primary ion energies appear quite difi'erently, the 
reconstructed distributions agree well with each other. The depth resolution in the reconstructed 
data is increased significantly and segregation of boron at the near surface side of the (5-layers is 
clearly shown. 
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The improvement of the depth resolution achieved in 
sputter profiling in general and secondary ion mass spec- 
trometry (SIMS) in particular has continued over the 
last ten years despite occasional predictions about the 
limit had been reached [Q. Nevertheless, the resolution 
achieved directly may not be adequate for future gener- 
ations of semiconductor device material, even if ultralow 
energies |^,|| or large cluster ions |^ are employed in the 
primary beam. 

The SIMS mass transport effects due to energy depo- 
sition and probe beam incorporation from the primary 
beam are well known — the measured profile is broaden 
and shifted from its true position . Although the use 
of low beam energies or high mass clusters can greatly al- 
leviate the effects, such mass transport effect is intrinsic 
to the SIMS depth profiling process and can still be ob- 
served even when the primary ion beam energy is as low 
as 250eV g or when using SF+ ions at 600eV @. The 
true spatial distributions remain to be reconstructed es- 
pecially when there is an abrupt interface or 6- doping 
present. 

In the concentration range of common practical inter- 
est in SIMS, e.g. the dilute limit where dopant concen- 
tration < 1%, there is a strictly proportional relation 
between the signal and the instantaneous surface concen- 
tration as well as between the primary ion flux density 
and erosion rate outside the pre-equilibrium region for 
a single matrix Where the essential physics of the 
analytical process is linear, deconvolution is the math- 
ematically correct method to recover and quantify the 
depth profiles |10|. The ideal SIMS signal at the depth 
z, Yq{z), can be expressed as a convolution of the true 
concentration distribution, C(z), with the SIMS instru- 
mental response function, R{z), 



Yoiz) = / C{z')R{z- z')dz' 



(1) 



since primary ion flux must be also proportional to the 
elapsed time (otherwise depth calibration becomes prob- 



lematic). Here we define Yo{z) and C{z) with the same 
dimension of concentration per unit length and R{z) nor- 
malised over the depth to simplify the equation and en- 
sure that sample mass is conserved. A slightly more so- 
phisticated model Yo{z) = J C{z')R{z — z\z)dz' might 
be used if the depth resolution was depth dependent. 

It would be a straightforward inverse problem to de- 
termine the true spatial distribution C{z) if the corre- 
sponding ideal SIMS signal Yq{z) could be measured. In 
fact, the measured SIMS signal, Y{z), is as usual a com- 
bination of the ideal signal and associated non-negligible 
noise component, Yn(z), 



Yiz)^Yo{z) + YNiz) 



(2) 



There is no obvious way to find C(z) from Y{z). Various 
methods have been used to obtain C{z) Q. Yet, the lack 
of objective evidence makes it very difficult to reconstruct 
the real features and separate them from the SIMS effect 
justifiably, e.g., to distinguish the segregation and diffu- 
sion occurring during growth at the interface of two dif- 
ferent materials from the SIMS atomic mixing. Moreover, 
the peculiar character of the SIMS depth profile data ev- 
erywhere positive, and large dynamic range (may span 10 
orders of magnitude overall and 4-6 orders for a particular 
species), requires a very careful and unbiased treatment. 
Manipulating data arbitrarily, such as placing a subjec- 
tive penalty on each change in slope or simply filtering 
certain range of frequency components, could seriously 
distort the final results and/or easily lead to unphysi- 
cal negative values. We believe that only the features 
with statistical evidence in the original data should be ex- 
tracted and a empirical deconvolution method | based 
on the maximum entropy (MaxEnt) principle ||l 21,131 ful- 
fils such a requirement. 

The success of a MaxEnt deconvolution method relies 
on the finding of a sample independent response function 
and a suitable noise model. Neither of them can be ob- 
tained for the SIMS process with sufficient accuracy from 
first principle calculations because of incomplete knowl- 
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edge. Our recent investigation shows that the noise in 
the SIMS depth profihng follows Poissonian statistics uni- 
versally [0 and the relation between the mean counts, 
Sm{z), and its corresponding standard deviation, (7s{z), 
is as{z) = Sm(z)i/2^ 

The response function should contain only the SIMS 
related information, i.e. the broadening, the shift and 
possibly the ion yield, but must not contain sample de- 
pendent features. Ideally, the response function is the 
transient measured from an infinitesimally thin layer, of- 
ten known as a (5-layer. However, such a S- layer is an 
abstraction and even if it were not, we have no means to 
recognise its existence, other than very locally. Moreover, 
as the intrinsic resolution of SIMS is improved by using 
sub-keV probes it is readily apparent that the real 
approximations to such structures which can be grown 
leave a measurable sample- related shape content in the 
SIMS profile due to statistical placement of atoms, seg- 
regation, and diffusion at the growth temperature, etc. 
Simply using such data as response function, deconvolu- 
tion will suppress real features in the depth profile, and 
produce artificial concentration slopes and unrcalistically 
small feature widths. 

In the following, we outline a method for extracting a 
sample independent SIMS response function for the case 
of boron in silicon sampled by oxygen beam for various 
primary beam energies. Subsequently, we demonstrate a 
MaxEnt deconvolution to reconstruct the dopant distri- 
bution from SIMS depth profiles using the corresponding 
response function and noise model. The results are then 
discussed. 

Experimental and theoretical studies explored the 
mass transport in the SIMS process [p^-p^ . It has been 
shown that the normalised response function can be rep- 
resented by the following form in several orders of mag- 
nitude [Rl: 



R{z) 



2(Ag + Ad) 

{[l-erf(Cg)]exp [z/Xg 
+ [1 + erf (Crf)] exp [-z/X^ + (a/A<j) 2/2] } 

where £,g = (z/a + a/ Xg)^^^ and = {zja ~ a jX^'^l'^ , a 
is the primitive standard deviation, Xg and A^ the growth 
and decay lengths. The smaller the Xg and Xd are, the 
sharper the Riz) will be. When they approach zero, the 
i?(2:) will degenerate to a Gaussian distribution with the 
deviation a . All these parameters of the response func- 
tion apparently depend on primary ion beam energy, iJp, 
and should be monotonically increasing functions of it. 
If we could find a perfect (5-layer, we might be able to 
fit the measured data with R(^z) to obtain these parame- 
ters for a certain Ep and use them directly to deconvolve 
other SIMS profiles. Therefore, the difficulty here is how 
to determine the cr, Xg and Ad from the measured data 



justifiably, since there are no other techniques which we 
can use to check the results. We have to substantiate our 
choices through statistical and trend analysis. 

Before we study the Bi^z) of a specimen in a crystalline 
substrate, we first look at the Riz) of the same spec- 
imen in a corresponding amorphous substrate. This is 
because that the amorphous (5-layers are grown at very 
low temperature where segregation is negligible and only 
broadening due to diffusion is present. Hence the amor- 
phous i?(z) should be qualitatively the same as the true 
crystalline one. We measured an amorphous boron b- 
layer in amorphous silicon grown at room temperature 
by molecular beam epitaxy (MBE) with the range of Ep 
from 335eV to llkeV and fitted the measured profiles 
with the R[z) in Eq. (3). The primary ions used are 
normally oxygen incident. The obtained parameters are 
plotted against E^ in Fig. 1. Within the error of mea- 
surements, we have cr""P'' = 0.86 -I- 0.27£;°-^^ A;;'"p'' = 
and y^J^'P^'- = 1.60i?p'^'^, where the lengths are in nm and 
Ep in keV. This reveals that the SIMS mass transport 
effect itself has no contribution to A^™^'* . We believe the 
same is true in crystalline case. Moreover, as the SIMS 
mass transport effect will be minimised at zero beam en- 
ergy, it is reasonable to think that the X^^p^ [g only due 
to the SIMS effect and the residual cr°'"P'' at Ep = 
comes entirely from the boron amorphous S- layer itself. 
This is consistent with the well known fact that a con- 
served diffusive point source has a Gaussian spatial dis- 
tribution. Consequently, the SIMS contribution to the 
amorphous primitive deviation, (j°™p'»^ can be obtained 
from cr™P'*(£;p)2 = a'''^P'\Epf - (t™p''(0)2. 

Crystalline R{z) is built up by fitting the SIMS data of 
a MBE grown boron (5-layer in a crystalline silicon sub- 
strate with the range of Ep from 250eV to llkeV under 
the same experimental condition as for the amorphous 
study. Considering that the bonding energy of atoms 
in a crystalline material is usually larger than its amor- 
phous counterpart, we expect that both a'^'^^"* and X'^J^'^^ 
are smaller than the amorphous ones and A^''*"'* should 
remain zero. The boron i5-layer structure grown by MBE 
normally has a segregated interface at the near surface 
side and a almost ideally abrupt interface at the other 
side. This will enable us to get a reliable X'^J^'^* from 
the fitting. For A^'"^'** , the fit shows no significant energy 
dependence and we therefore take it as intrinsic to the 
sample. The results are also shown in Fig. 1. We obtain 
that a^'-y* = 0.27 + O.SdE^ '^^, and A^^'* = 1.39E°-^'^. 
These parameters are indeed smaller than the amorphous 
ones, and again X^jy^^ vanishes as Ep approaches zero. 
The (j'^^y^* can be partly affected by the sample structure 
in use and, indeed, shows a finite intercept. Therefore, 
we take the intercept away and calculate the SIMS part, 
cr'=''2'^*, as afy'\EpY = a^'^v'^Epf ~ a^'^v'^Qf. 

Using the obtained parameters, we are now able to cal- 
culate sample independent SIMS response functions for 
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various Ep and use them with our noise model to de- 
convolve some measured SIMS depth profiles by MaxEnt 
method. A multiple boron 5-layer structure in silicon was 
grown at a constant 450C by MBE with the arrangement 
of, in turn from the surface, a pair of (5-layers 2nm apart, 
then another pair of J-layers 5nm apart, etc. The in- 
tended boron concentration was 1 x lO^Oatoms/cm"^ for 
all the layers, which is well above the bulk solid solubility 
limit ijl^,^ . We profiled the sample using a normally in- 
cident oxygen beam at Ep — 0.25, 0.50, 1.0, 2.0, 4.0, and 
6.0keV, respectively. The depth calibration of profiles 
in silicon using oxygen primary ions has been discussed 
and clarified recently ||2l| . In our case, since the part 
we study is outside the pre-equilibrium region and much 
longer than the transition depth we can align centroids 
of the profiles based on the principle that distance be- 
tween centroids will not be changed by any convolution 
process. The MaxEnt algorithm used in our calculation 
is from Ref. |Q. The measured profiles after conven- 
tional calibration and corresponding deconvolved spatial 
distributions for the first two pairs of boron (5-layers are 
shown in Fig. 2(a) and (b), respectively. 



structed spatial distributions agree well with each other. 
The area under each feature is the same before and after 
deconvolution, i.e., the number of the boron atoms for 
each feature is conserved. It is obvious that the depth 
resolution of the reconstructed spatial distributions has 
increased significantly. The J-layers at z — 87nm which 
were separated by 2nm can just be distinguished at 
Ep = 2keV while the 5nm pair at z — 135nm are well 
separated at the maximum Ep = 6keV we used. Note 
that the deconvolved features become smoother as Ep 
increases, i.e., better depth resolution can be achieved 
as Ep gets lower. From the reconstructed dopant dis- 
tribution in Fig. 2(b), it is clearly shown that there is 
considerable boron segregation on the near surface side 
of the layers. This kind of feature would be automat- 
ically eliminated if a simplistic R{z) were used. There 
are some unexpected small concentration spikes in the 
reconstructed data at level of two to three orders lower 
than the peak height. We have no objective criteria to 
confirm their existence or take them away so they remain 
whether one likes them or not. 
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FIG. 1. The response function parameters, a, Xg, and Ad 
vs. primary beam energy Ep. The open and solid symbols 
as well as the dashed and dotted lines are for the amorphous 
and crystalline parameters, respectively. The lines are fitted 
with form a + b Ep for each of the parameters. 



From Fig. 2(a) and (b), we can see that although the 
original SIMS profiles appear quite differently, the recon- 
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FIG. 2. (a) The SIMS depth profiles at different primary 
beam energy Ep for the boron 5-layers in crystalline silicon 
substrate and (b) the corresponding spatial distributions re- 
constructed with the SIMS response function through a Max- 
Ent deconvolution. 



Using Eqs. (I) and (2), we can easily work out the cor- 
responding generalised Rayleigh limit of depth resolution, 
Az, for two adjacent ideal 5-layers after deconvolution if 
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we assume that Az is limited only by the SIMS noise: 

^z{Ep) » ^UE^) (4) 

where we take \g = 0, 3> (Ts, and Yc is the origi- 
nal measured counts. Estimating with the typical peak 
counts Yc ~ 10'* in our experiment and the from 
Fig. 1, we find the resolution we have achieved in the 
above deconvolution is within the limit. 

Maximising the entropy of spatial distribution gives 
us the most likely deconvolved solution. It has a global 
tendency to spread the solution onto the whole space 
range within a given noise deviation when total concen- 
tration is fixed. However, there is no constraint on the 
local changes of the distribution. For example, two spikes 
in the distribution contribute the same entropy no mat- 
ter whether they are next to each other or not. This, 
unlike local constrains such as limiting the derivative of 
spatial distribution, makes possible the reconstruction of 
some abrupt features, e.g. <5-layers, superlattice struc- 
tures, step doping materials, etc. 
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FIG. 3. The Fourier power spectrum, of the mea- 

sured SIMS signal Y[z) vs. wave number, k. The data are 
calculated from the 0.25keV SIMS boron profile and the con- 
volved one from the corresponding reconstructed distribution. 

To have some further understanding on the MaxEnt 
deconvolution method, we compare the frequency compo- 
nents of the SIMS measured profile and the correspond- 
ing convolved one calculated from the reconstructed dis- 
tribution. Fig. 3 shows the power spectra for the Ep ~ 
0.25keV case. Although the MaxEnt method only ma- 
nipulates data in real space, the background noise fre- 
quency components in the convolved profile are clearly 



suppressed by 6-7 orders. It is clearly shown that the 
high frequency features in the range of 0.75 to 1.50nm~* 
are retained rather than eliminated if a conventional 1// 
noise subtraction method were used. Profiles of other Ep 
have the similar results. This is consistent with Shan- 
non's entropy loss theorem for signal passed through lin- 
ear filters Consequently, the MaxEnt method is not 
only able to reduce the background noise level without 
using any artificial windows in real space or filters in fre- 
quency space but also capable of retaining sharp features 
which are statistically significant. 

We should emphasis that although the MaxEnt 
method does improve the depth resolution and recovery 
sharp features, it only provides us with the statistically 
evident information which the original SIMS profile con- 
tains. The MaxEnt deconvolution is not an alternative to 
instrumental improvement, such as achieving lower pri- 
mary beam energy and higher corresponding beam cur- 
rent. For example, one cannot distinguish whether there 
are two close 5-layers or a single 5- layer of higher concen- 
tration if the profile is measured with very high beam en- 
ergies. This is because high energy profiling will not only 
have large A^ but also lead to large depth increment and 
low sampling density, which will limit the information in 
the data at the first place. To reduce the depth increment 
for higher sampling density in this situation requires sig- 
nificant reductions in the primary beam current because 
of the high sputter yield at high beam energy. This will 
lead to lower values for Yc. Therefore, the signal- noise 
ratio will decrease, which limits the potential increase of 
depth resolution from the deconvolution. 

We describe a procedure to obtain sample independent 
response function which is suitable for use where the in- 
trinsic SIMS effect is so small that sample related features 
can be seen. Using this together with an empirically 
determined noise model, MaxEnt deconvolution is per- 
formed to reconstruct self-consistent dopant distributions 
from the SIMS depth profiles obtained at different beam 
energies. There are no adjustable parameters involved 
in obtaining the response function and deconvolving the 
profiles. The depth resolution of the reconstructed dis- 
tributions has been greatly improved and the segregation 
on the near surface sides is clearly demonstrated. With 
reconstructed true spatial distributions, further quanti- 
tative investigations on interface segregation, atomic dif- 
fusion, and related problems will be straightforward. Our 
method can be used to study other dopants, and the Max- 
Ent formalism may be extended to models other than 
convolution. 
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